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TEE ELLIPTIC tflKG BASED OH THE POTENTIAL THEORY* 
By Zlaus Krienes 



SUMMARY 



The present report deals with the elliptical wing in 
straight and angular flow on the hasis of the potential 
theory. Conformahly to the theory of first approximation 
upon which the calculation rests, the |^nown requirements 
regarding the shape of the surface and its angle of attack 
must be met. A further condition is that the slope of the 
surface toward the streamlines must be a continuously dif- 
ferentiable function of the points of the surface. If 
this is not the case, in a given example (for instance, 
by aileron deflection or wing dihedral - the latter being 
of importance in sideslips), the discontinuities must be 
replaced by suitable rounding off. In general, the cal- 
culation of a given elliptic surface requires a series of 
infinitely many potential functions, the coefficients of 
which are afforded from linear infinite systems of equa- 
tions. The expansion is stopped with a certain term, de- 
pending upon the degree of accuracy desired. Its effect 
on the integral quantities, lift and lift moment, is prac- 
tically negligible. An immediate prediction of the in- 
duced drag is ruled out, since it would involve all the 
coefficients of the infinite number of potential functions^ 
Otherwise, the lift distribution at the wing tips does not 
approach zero or the downwash becomes infinite, which is 
due to the fact that the load distribution of the lifting 
line is developed here by spherical functions (equation 
(80)) which do not approach zero at the wing tips as do 
the trigonometric functions employed elsewhere. On the 
wing in sideslip, which can be summarily replaced by a 
lifting liiie, the so-called parasite drag (reference 2) 

= - / /(pu - Pob) Y dy 



*"Die eiliptische Tragflache auf potentialtheore t ischer 
Grundlage." Z.f.a.M.M., vol. 20, no. 2, April 1940, 
pp* 65-88. 
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would have to "be defined fir s t and the suction force on 
the leading edge suhtracted therefrom, where, however, 
extrapolations are recommended hecause of the finite num- 
ber of computed coefficients. Even the resulting pres- 
sure distrihution is only conditionally valid "by few ex- 
pansion terms, especially near the wing tips. 

It msiy he mentioned that the computed potential and 

doxfnwash functions change on transition to > 0 

into Xinner*s functions for the circular wing. 

A large portion of the computations were made on the 
calculating machine, the accuracy of the slide rule heing 
insufficient in che calculation of the elliptic integral 
for higher n. 



IHIEODUCTION 



This article is intended as a contrihution to the 
theory of the lifting surface. The aerodynamics of the 
elliptic wing in straight and ohlique flow are explored 
on the hasis of the potential theory. The foundation of 
the calculation is the linearized theory of the acceler- 
ation potential (references 1 and 2) in which all small' 
quantities of higher order are disregarded. This affords 
the follo\ving simplifications: 

1. The z coordinate of every wing point is neglected, 

i.e,, the variation of the potential function 
corresponding to the pressure pump on the sur- 
face is situated on the base ellipsoid (fig. l) . 

2. The streamlines, along which the convective inte- 

gration of the acceleration is effected, are 
straight lines parallel to the direction of the 
stream* 

In the cas e of the elliptic boundary, solutions of 
Laplace^s differential equation A 0, as products of 
Lame's function, are. known. 

The acceleration TDOtential and Lame^s f tine tions ♦ - 
Suppose that the elliptic \ving is in a stationary parallel 
flow with the velocity V. She fluid is homogeneous , in- 
compressible , f ric t lonless, and * not "sub jected tp gravity 
and non-vortical outside the lifting surface and the shed- 
ding vortex band. Then: 
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— = - _ _ grad p 
■ dt p 

The pressure p is expressed the potential func- 
tion ^ 

P - Poa = - P V« ^ (2) 

where Po, is the static pressure at infinity. Junction 
\i/ satisfies Laplace's differential eq.uation 



At}/ = 



^^^^^.0 (3) 



-V 2 -> 2 -.2 

oy, 0 2 



and may be visualized as being the result of a superpo- 
sition of sources and sinks of intensity cr(xj,, yj, ♦ Zj) 

on every point of the surface. 

^!/(x.y,2) - /7'ar(xy,yj,,2y) ^/^|)dP (4) 

/../ 

where n denotes the direction of the normals of the sur- 
face in {r.jf7^f^j) and 



the distance of the starting point (x,y,2) from (^c^, 
y-,2_). The integral can "be exchanged for one taken over 
the surface of the ground ellipse. 

\i/(x.y,z) = /)'^cj(xp,yp)^f i^) dxj, dy^ (4a) 

The mathematical treatment of the present case he- 
comes possible by the introduction of ellipsoidal-hyper- 
boloidal coordinates (reference 3). The semi-axes of th e 
base ellipsoid being c in direction y and ca/1 - 
in the direction of x (so that 2ko becomes the distance 
of the aerodynamic centers on axis y) the new system of 
coordinates is: 

^ y_^_ . z? ^2 



3 2 
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The solution of the equations leads to: 

^ ^ c r^. n — ~ ^k""' „ „ .... 



Z s 



, ° - 1 ^1 - yi - V' 



f (6) 



The surfaces p, jj. , and v = const., respectively, 
are confocal surfaces of the second degree; P = 1 yields 
the elliptic base surface, p, = 1 the plane 2 = 0 out- 
side of the elliptic disk. The surface element of the 
base ellipsoid is 

dxdy^c^ >-- -v dudv (7) 



besides which 



^1 - IJ,^ -s/l - =^ 



y- - 



(8) 



^/l - 1 - ^'^ 

is valid for p = 1. 

The introduction of new coordinates u, v, w 

V(u) = p^ - id + K^). Y(v) = - i(l + K^). 

Y(w) = - i(l +k2) (9) 

o 

for p, H-* i; means of v/eiers trass ' V function (ref- 

erence 4) , appearing ia equation (lO) 



^11:^1 = 2 V(^(u) . ei)(7(u) - e2)(7(u) - 63) (10) 
d u 

where the quantities 

ei = -^(2 -k"") ; 63 = ^(SK^^ . 1) ; e 3 ^ i( 1 + K^) ; 

ei + 63 + 63 = 0 (11) 



i . e • , 
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63 « 1 



are posted, gives 



u = / 



d p 



d p 



(1-2) 



'"4 S/y/p - eiVp - ea^p - 63 

00 

and Laplace's equation reads 

[V(v) - v(w)]|!| t [^(w) - Y(u)] |!|+[ Y(u) .^ycv)]!!!', 0 (13) 



6v' 



Posting the solution in the form 

\i;(u,v,w) = 3(u)E(T)il(w) 



(14) 



gives, for each of the three functions 2, Laplace's differ- 
ential equation 



d^S(u) 
d u'^ 



= [A + B Y(u) ] E(u) 



(15) 



with the separation constants A and B. After posting 



B = n(n + 1) and ■'"n = - 



1 + K 



8 3 



and^again introducing M-t V, p "by means of equation (9), 
Lame's differential equation (reference 3. yol. I, p. 359) 

reads 

(^.*-K^)(^^.l)i!^ii^+».(s^.« - K« - i)Hn!kl ^ 

dn^ d|A > 

+ [(l+K^>„~n(n+l)M.^j3a''(M.) '0 J (15a) 

For 3 = n(n + l) there are precisel:;- two 2tt + 1 
values Vn, for which E^^^dj.) has the form ( reference 3 , p, 360) 



(16) 
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On^(|i) aoM.^""^ i""^2-^3 + . . is an even polynomial in p. 
of degree n - - - €3. The values are dissimi- 

lar solutions of an algelDraic equation, resulting from 
the condition that On (m*) is a polynomial. V/ith a view 
to calculations later on, En®(v) is defined as follov/s; 

i K. ] Ik J ^ 
= A(p^i(-cos (p)^(ain cp)^=^0a°(q5) (16a) 

■by putting 

- = sin (p ; 'Z'^" ~ '-^"^ = -cos ; A^> = -/l - sin^cp (17) 
The solut ion 

M/'n"'(**.V.p) E/(^)En'"{v)En'"(p) (18) 

achieved x^rith equation (13) can also te represented in 
different form (reference 5), Tor, on denoting the zero 
places of the polynomial O^^vM') with Pg^t is readily 

seen from equations (5) and (6) that 

s^'^(^)s/(v)s/(p). 

One of the factors contained in the parentheses 
is selected and the product is formed over all zero 

places "beonging to Ov^ (m») . The zero places pg can 
also oe determined direct hy applying Laplace's operator 

7^^ Ti^ / \ 

— p- + 4- to the right side of ecuation (19) and 
3x*^ oy^ dz*^ 

making the result equal zero. The system of equations for 
is then as follov-s: 



3^2 3^3 3^1 4 8=1, 2, ..•^ 

"3 5-=0 (20) 



Ps Ps Ps " Pq m= n- €, - e^- € 



KArCA Technical Memorandum No. 971 



7 



The potential functions in the respective form of 
(18) and (l9) are so-called "inner" solutions of the po- 
tential equation, not suitable for our purposes, since we 
require potential functions which ordinarily disappear at 
infinity. These are secured hy taking the Lame function 
of the second type in the variable p. This is the solu- 
tion of Lame's equation (l5a)» which, for p — > oo as 

^^^J^ —^0. It has the form 

pn+ 1 - 

00 

^V(P) = s/(p) f (21) 

J [Ea Cp)J ^/p2 - 1 Vp2 - k3 
P 

The integral can always be reduced to elliptic in- 
tegrals of the first and second categories only. The 
aspect of the outer solution of the potential equation 
is then as follows: 



= E/(^x)E^'»(v)s/(p) T— — 



(22) 



P 

and eq,uation (19) yields 



, m, V const J , ^ I TT _^ 

° I z zx ) Ps - Ps 




s 

d o 



(23) 



From the representation of the potential function by 
(4a) as source-sink superposition on the elliptic disk, it 
is appc.rent that the potential functions in the plane 
z = 0 in the outside zone of the disk must be zero. Hence 

*4^a^ containing factor z must be taken according to (23) 
i.e., only such Lame functions as are of the form (cf. (6) 
and (16)): 



E/(^^) = ^A - M.^ Mj^^(i^) (24) 

Lam^^s functions have orthogonality characteristics 
similar to those of the spherical functions (references, 
vol, I, pp. 369 and 379): 
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^ A ^^ ■ r Ot if m 1^ t . 

Jhc - V® Vr^ - , if m = t 

•../ V 

r 0, If n ?^ s 

^ ■ °' * (26) 

1^, if n = s 

and m = t 

i'l-pi connection with etiuation (4a) is established "by- 
having recQUl -ae to the following represen tation of 

^ (a =y(x - xp)2 + (y - yj)2 + (z 2j,)2) (reference 3, v. I.I, 

p. 172) : 

n=o in = 1 

(27) 

whence 



3 2p \ R/ 2j,=0 



leaving the summation over 3j^°^(p,) in the form (24) to 
he effected. Then hy assuming that the source-sink distri- 
"bution on the elliptic disk approaches zero o n the e dge with 
the root from, the edge distailce , i-e., with 1 - p.^,^ , 

0"(xp,y-j) can he developed cdnformahly to products of Lame^s 
:function of the type (24) i 

CO 

rr(Xj,,yy) = As* .(^F)^s*(^r5 <28) 

after which the.' Toi'mulat ion of the. integral (4a) gives, 
hased on the orthogonality of Lame's products (cf. (26)), 
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^H\^.V.p)^^/l-K^ 2 A/l/Mj,^(l)E/(^.)E^"^(a;)Ja^(p) (29) 

Sq^uation (29) proves the potential function (4a) to 
"be a certain linear combination of the functions ^j^^ t 
which are analyzed next. 

Hepresentation of the potential function ^^^^ as a 
definite integral .- If 



3 r- - ■ 

X = 2 i ^ 'A(t)-ea W(u)-ea 'A(v)-ea v'v(w)-. 



^^a-3)(e,.e, ^^^^ 

(e^, defined \>y equation (H)) or if (6) and (9) are taken 
into account 



X = ^''(O - 62 V^(t) -63 V^(t) ~e, Z 

Ve^-elve^-es^ ° */ ©3 -©iV 63 -eg ® '/©i -eg*/ ex -63 



it can "be shown that 



\|;j,'"(x.y,2) = / <in(X) En'°(t) dt (32) 

t-'TTr 



i ^/ 63-63 i/eg-es 

X polynomial [Y(t)] ^a^lJ ^^S) 



if a loop about the points of the complex t plane cor- 
responding to X = il is taken as integration path from 

toward TTg. Qn^^^ satisfies Legendre's differential 
eauation 



(1-X«) il^ia.2X ^(^+1) n (X) = 0 (34) 

dX dX 

Ej^^(t) complies with Lame^s differential equation 
(15). OJhen 
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. a^^XX)^ ^ ^ l')[Y(u) - V(t)] q^(X) (35) 

for any two variables each of t, U| v, w; i»e., 
S/(t)[^^^-{^ia + ^(n+l)Y(u)}ci^(X)] 

= 3 »(t) . 0 (X) £la!lil (36) 



and 



3u 



2 V(u,v.w)-. [Ao + n(n+ l)Y(u)3 '-It (u.v.w) 



1 r."(o^^-4.(x)i.Ia!lI>r' (3,) 



2TTi L ^ at n d t 



-J IT, 



The same holds true if u is replaced by v or w. 
The integration path is next so chosen that 



i E^n^Ct) MaiAl- q^(X) ^- V^*^ 

L at dt 



"2 

at dt J TT 1 



i.e., •(i/jj^ in every one of the three variables u, v, w, 
satisfies Lame's equation (15) and is accordingly a third 
representation of the potential function defined by (23). 
The points X = ±1 in plane t are given by 

t=v+w+u (X= -1) ; t=v+w-u (X= +1) (38) 

On the elliptic disk (p = l) , •V(u) s ei; i.e., 
1 

u= r ^ P =~u> (39) 

J^Jp^ - 1 Vp^ - 

The potential function ^^^^^ is cited as. an example.. 
The sole Lame function of t he fir st type and first degree 
equipped with the factor #/p® - 1 is: 



3i^(p) =yp2 - 1 
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According to. (22) and (33), respectively, the poten- 
tial function then reads, respectively J 

dp 



/— — * P 

^!/ 1 (p.|Jk,v)=<v/p^-l'/l-»AVl-v^ / — — y-. r 



and 



Jl-K^ p dp 



' C / (p2. 1) Vp^- 1 VP - tt^ 



Lift and lift moaents. - The lift is given hy 
A = / / (p - Po.h)'i3cdy = pvV/(^koh -^l'u)<i3cdy (40) 



I ell 
wherehy 



n ip=i ° " ^ Mj^"^(i)yi^ 



■4/ 



Based on the or the goaality of Lame's functions, only 
. contril>utes to the total lift 



1 



A = 1 £ 7^ F (42) 
3 2 ell 

furnishes the pitching moment about the y a.xis 

M = ^ c ^1 ^ I 7^ Fell (43) 
^ rolliixg moment about the x axis 

= ^ ° 2 ^ell (44) 
(The negative index refers to the odd functions in y,) 

The eliiDtic wing in straight flow* - Assume the el~ 
liptic wing in a stream in direction of the positive x 
axis with velocity V, Kow equation. (l) enables tTie cal- 
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culation of the velocities induced the pressure poten- 
tial in space and especially on the lifting surface. 
The z component of equation (1) reads in the stationary 
case 

8x dy 3z p Sz oz 

Small quantities of higher order are disregarded, 

i • e • » 

^ = y2±. (46) 
Ox dz 

The z component w of the velocity vector w is 
hereafter called "downwash" for short. The downwash on 
the elliptic surface is obtained by integration of equa- 
tion of equation (46) for z = 0 and y < c over x: 

X 

w ^ r M d X (47) 

V / dz 



The calculation of the integral is readily secured "by hav- 
ing recourse to the representation (32) of the potential 

function ^' . After formulating hy differentia- ■ 

n oz 

tion "below the integral sign, equation (47) gives 



^ r r JL- / Vv(t) - e2 X 

2Tri / / oz^^\ r — r — " c 

J J y Ves - e^^ V^s - ^3 



^^' ^ ^ - 1 + ; ^ En^t)dtdx 

- e^Ves - 62 ° ft/ei - eg Ve^ - 63 (48)J 



32 ~ dX 8z - dX Ve, - e^v^e, - ° 



5 ^n(X) . d^nCX) SX _ iCjnW -Jyjt) - es 1 

Sz dX dx dX Ves - ei ^/e^ - 63 
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that is I 

(49) 

^2 3 X yY(t) - eg 

which, when inserted in equation (48) and followed hy in- 
tegration over X, affords since lim ftjj(X) = 0 for 
X — > ^ - 



'•^(X) is given Tsy IJeumann's representation. 

+1 

1 r p (Y) 
"^-1 .... 

Similarly X(t) is expressed with 



ve3~eiV®a~63 Ves-CiVes-eg vei-egvei-eg 
whence (50) becoDes 



Ik / ^ / — ° — ~======iSSn (t)dsdt (53) 

V 2TTi / 2 / X(t) -I(s) ds VV(t)-e3 

Q 

The integrand has poles at t - s, because X ( t) 
Y(s) = 0 and at t = iv^ + ^2* where V7( t) - e^ has a 
simple zero place. The "behavior of the denominator near 
the zero place is defined hy Taylor expansion. It affords 

y;Ur^3 = Ct - (a,, (^^/P(t)-ea)^^^^_^^^,^ + .. 



- - [ t - (lOi + i 2)] ^/eg ^ e^ '/es -03 + • • . (54) 
taking into accouKt equation (l), as well as 

X(t) - Y(s) = (t- s)^ (t - s) ... (55) 

V 8t/t=s °* 
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according to equation (52). Then the integration of 
ssv+w + Ti as far as s =.v + w; -u corresponding to 
Y = -1 to T = +1 (equation (38)):gives 



1 n Pn(I(s)) iI.ft^(xK^l >a;3)) (56) 

2 / X(»^i + i tU3)-y(s) ds ^ 



By 

Mb) - 



/ P^(y(s)) , ^'^"^^ e/(s) d s (57) 
J yY(s)-e2 

it is to be noted that the integrand for p = 1. i.e., 
u = -ujj^ has the period 2 a'j, so that the integration 
path can be shifted until s proceeds from i ujg - uj^, 
to i a)2 + (fig. 2), which corresponds to ^ " ~ ^ 
to when ^^(s) - e, = - K sin e , that is, ds = 



||. (Ae = Vl - K 3 giaS e). Moreover, let 

y 



CVl - K' 

80 that, because of 



= I. - = (58) 



X(a)i + iujs) = ==^1—:^== J = n (^(101+ iujg) s eg) (59) 

and, according to (33): 

Ea''(Wi + iWa) i Bn^(X)^^^ , (60) 

the downwash function on the elliptic disk becomes 

• • ■ s ■• 

^= Sn"'(K)qj,(,l)-i /^'PnUcose + ^3ixic)M„^(c)A€-^ (61) 

V I ■ .'" . '. CO S€ 



2 



The coefficients ohtained in. the polynomial' of i 
and Ti are complete elliptic integrals of the first and 
second categories. The calculation of the downwash function 
in the case of n = 1 is expressed as: 
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Pi(T) = T, ^i(X) = I m 1^ - 1 



From (61) follows 



''^''^^ '"^^ = Jl - k2 (<i(t») - i /* (I COS€ + Tl Sin € d€ 
V / C OS € 



that is , 



g he lifting surface .- The shape of the surface is 

given "by 

z = 2(x,y) 

The slope of the surface in x direction must agree with 
the direction of the flow at the same point, that is, 

Q ^^^>y) = >^(^fy), (62) 
• ox V 

from which follows, for z = 2(x,y) 

X 

z(x,y) = i y^w(x,y)dx (63) 

the lower integration limit being arbitrary; we equate.lt 
to zero and add to the value of the integral an arbitrary 
function in y: 

X 

2(xty) ^ P w(x,y)dx + g(y) (63a) 

/ 

• 0 

The lifting line, the induced drag, and the suction 
f orce>- W'e merely refer to the corresponding chapters of 
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Kinner's report (reference 2), where these problems are 
treated in detail. The results are readily applicable 
to the elliptic disk. The lifting line, by x/hich the 
lifting-surface is assumed replaced ,• is obtained by cor- 
relating the lift elements through integration parallel 
to the X axis: 

|^(Pu-Pob)ix = 2p V= r^n^'U,^ 



a^"'(Ti) 



E 



It affords, :for Instance, 



aii(n) = 4Tr ^ 7* c -Jl - 



?fi(Tl) - P8(r.) 



agi(Ti) s 0 



> 



a4^('n) s 0 



P2(T1) - P4(ti) 



(64) 



The -potentiel function of the second type ^-* The fore- 
going potential functions lend themselves in any way to 
linear combination and yield the corresponding linear com- 
binations for lift, lift moments, and downvash w. The po- 
tential functions dealt with so far afford the aerodynamic 
q^uantities of a correspondingly curved wing hy shock-free 
entry of flow, that is, at a certain angle of attack where 
no flow around the leading edge occurs. The arbitrary 
angle of attack is obtained by superposing a flat elliptic 
disk with its flo\ir, where, as is known, the leading edge 
is suction edge, that is, the lift density approaches in- 
finity. All the potential functions of the first type 
approach zero, however, on the disk edge with the root 
from the edge distance; hence the task of finding poten- 
tial functions that have these qualities. They are 
achieved by applying on the potential functions of the 
first type at constant x, y, z, and. the following 

boundary transition (reference 2): 



(66) 
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These potential functions- possess the quality 

of hecoming inf inite on the whole horder of the disk* 
Because of the condition of smooth efflux on the trailing 
edge, it later is necessary to comhine the functions 
linearly. 

The downwash function of the second type ** In con- 
formity with equation ( 65) , the dovynwash functions ate 

obtained hy applying the operator — ^ l^c^to equation 

^n- X do 

(61), while oDserving the interrelationship: 

The same applies to T^^i'^) aad <^(t[). Equation 
(61) then gives +^ 

V ^ d ti 2 / 4 y a cos € 

3 (66) 

In dealing with the .second fundamental problem, that 
is, in the calculation of a prescribed wing, potential 
functions are used, the downwas.h functions of which are 
independent of x on the disk, According to (46) , this 

• 3w 
implies, since — = 0, that 

8x 



s 0 (67) 



With coefficient 1?^° still to oe defined, we put 

0^(x.y,z) = £ 1,/ ^/(x,y.z) $.^(x,y,z)= Z^^'"'^^-'^U,7 ,z) 
m -•m 

(68) 

Correspondingly it is; 

= Z I,/ .^^11, w.,. _Z (68a) 
Wj^(x,y) is designated as downwash function of the poten^ 
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tial function of the secohd^iyjJ^^* <^n^» ^ (68a) 
and (65) , it is : 



(69) 



2 



The coefficients t^^^ are now so defined that w^^ 
is a function of only, The first, term in (69) already 

depends on n only; hence th^' second term itself, which 
usually depends on i also, must he a function of ti only. 

d p ^ y ) 

— is a polynomial in Y with terms of the form 

d Y 

Y ■t' = |_ ^cos e + n sm ej 
hence is a sum of terms of the form 

(cosc)''"^^^'' (sin c)^ |n-2P-a (^q) 

For the following arguments it is assumed that Mj^^(c) 
is even in € ; so that all terms with odd powers of sin c 

disappear in the integration from to But if a 

is even, (sin e) = 1 + sum of cos terms. V/hen this is 
written in (70), y^-sp consists of the following sum- 
mands : 

(cos €)^"^<1 ^n-2p.a (^q^) 
and the condition that 

n-2p-a<n-2q<n-2p (71) 

Putting 

cr = ^ and - ^ ^ respectively (a* is an 

^ : : .5 :: integer) (72) 

and stipulating that 

r (cose)^-^^ 2 t^%°(c) Ac_i_i_=0; q= (7S) 

'in """^ ^ 

~2 1 , 2 , . . . (cr - 1) 
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causes all terms containing the powers of cos €, i.e., 
powers of I , to disappear. The sole nondisappear ing soim- 

mand of Y^'^P is obtained when n - 2 4 = 0, that is, 
when first n is even and, a'ccording to (7l) n - 2p - a = 
this then reads, according to (70a), -n^^^P and depends 
no longer on so that it can he put hefore the integral, 

Comhining the summands hefore the integral, which now has 
the same value for every p, there is ohtained conformahly 

d p ( y ^ 

to in (69) for the integral 

d n 2 I m cos € 

"J- E 

(The integral disappears for odd n.) 

Mj^"^(€) contains the factor sin €; hence it is odd 
in €. Considerations corresponding to the foregoing then 
give the condition 

(cos c)^-^<i-' sin € L bn-"" Ma-^(e) Ae^-^ = 0^ 

^ " cos € 

q = 1.2 (T - 1); T= I and - ((74) 

respectively 

The value of the integral in (69) is other than zero 
only if n is odd. Summed up, it affords, "by attention 

to S|r(»^) = Ei?+i(ic) = 0 (equation (16)): 

i Wgr+iCri) = ksr+i — 5 -^v/^rVn) = igr 

V d T| V a T\ 

1 V/ \ . d P2r(Ti) 1 / \ 1 ^ 42r-x(^) 
i w.(2r+i)(Tl) - 02r+i — ' 7 w.8r(^) = ^sr " 

V a T\ V Ci x\ 

whereby 
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2 7 a 



? bir Mir (€) -ff 



cos c 




V (75) 



d € 



sine 2^ b;?+i «2r+i(OA€ ^^^ ^ 



1)5? Ei? (k) 



t^"* and 'bjj'^ satisfy equations (73) and (74), re- 
spectively. They are cj - 1 homogeneous equations for the 
c unknown "bn°. and T - 1 equations for the T unknown 
"bn"*"", respectively, which can te determined therefrom up 
to a common constant factor. The latter is so chosen that 



l=sr+i = lar = - Vl - 



•ar 



2r+x 



= 1 



(75a) 



whence 



(76) 



lw.(3r+i)(n)=^^^^; iw-3r(Ti)=.^^i-^5£Z2i!Ll 

Lift, lift moments, and thp lifting l ines of the ^- 
tential function of %iie second type .- These quantities are 

ohtained "by the application of the operator — - — c^ 

c^"^ dc 

to the corresponding quantities of the potential factor of 
the first type (42, 43, 44, 64). It is pointed out that, 
during the differentiation, the areal content of the 
ellipse Fell emhodies the factor c^. Then, "bear in 
mind that (equations (68), (75), (76a)) 



$ gives the lift A = 8 f P 



ell 



(77) 
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$g the pitching moment M» 



O Pa 

$-3 the rolling moment L = — c ^ V* 

3 ^ 



•J O 



sin^e d c 



(79) 



and for the lifting lines 



ai(T») = 4 ,TT |. c-Zl - k2 P (t^) 
aa(T\) s 0 



a3(Ti) = 4 TT I V* 6^1 - ic^ PgCn) 



a4(Tl) 5 0 



(80) 



The second fundamental problem of airfoil theory , - 
This involves the calculation of the aerodynamic quanti- 
ties of any given elliptic wing by means of the foregoing 
potential functions of the first and second types. The 
"boundary condition t'o be met is given by 

X 

(62a) 



which, differentiated, leads to (46): 



1^!^ = V 1^ (46) 

OX 02- . , - .. 

Posting ^|/ as a linear combination of ..potential 
functions of the first and second types, leaves 

= X a/ M/^° .+ E + 1 Da (SD 



leaves on the disk (p = l) 



8z 



= S E a^ 
n m n 



c z 



since 



O z 



p=l 



p=l 



• 0 (81a) 
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where then, however, 



p=i 



^A - M.2 -n/i - V3 Ld( -v/pS - 1)_ 



(82) 



p=x 



dw 



iiow, if -g^ is' developed according to the Lame 
products 



OX c n m " ^ " 



then, because of equations (46), (81a), and (82), based 
upon the orthogonality of Lame's product. 



d P 



m 



(83) 



This defines the coefficients ajj^ of the potential func- 
tion of the first type. On expanding 2 = z{ l,r\) in, a 
series of i and , the a^* ^s are determined hy a 
comparison of coefficients in all terms affected with pow- 
ers of i in thfe form of a screen meithod (reference 2). 

Prom equation (46) 



V n ^m 



m 



^dx= iw(r,) 
02 V 



(84) 



r-00 



is now only a function of tj* This residuary condition 

is complied with through a suitable linear combination of 

the downward function of the second type, v/hich depends 
BOlelyon t^. -The' condition reads 



^ Wn(Tl) + S Dn w.j,(r,) = w(n) 



(85) 



This equation is integrated from ti =^ 0 to tj , and 
then - multiplied by Pga- i('n) and PsYCr,) , respectively 
integrated again from ti » -1 to +1, and so yields 
through the intermediary o.f the int e„r-relations 



NACA Technical Hem or an dim No. 971 



23 



/ i*2r-i(n) P2a-i(^) d T) =S2(2a-l)+l 4a-l 
•-^1 . . Lo if r ?t a 



/ 



+1 

r> I 

'^ar(n) Psor-idl) dii = ; 

a^2a- 1) - r(2r + l) 

-1 



two infinite systems of equations for the coefficients 

+1 r\' 



and Dj^ 



/l- K 



2 



4a-l 7jJ 



r»o a(2a-l)-r(2r+l) ^a-l ^jj 

-10 +1 n. a =1,2,... 

Dsv+i — i— + S ^ ""^'^ " — \ r rw'(Ti')dT)«PaY(Tl)dn 

^^^4^+1 r=i ' Y(2Y+l)-r(2r-l) V / / ^ 



X86) 



The outflow condition** * A-part from the compliance of 
the equation systems (86), the coefficients Cj^ and Dj^ 
must also satisfy the outflow condition, that is , they 
must, he so chosen that the lift density disappears on the 
'trailing edge. The potential function of th^ first type 
satisfies this condition without that, since they disap- 
pear on the whole edge. The behavior of on the 
boundary is known with the application of the differentia- 
tion process (65) to (41), 



C^n"'^p=i= Mn'"^^) + ((y^)) 

M,a:i yc^ - ys - — 5^ 



(The term ( (Vl - M*^) ) disappears with the root from the 
edge distance*) We post, respectively, 

Z = M^(v) and g V^^^iSP) = M^Cqp) (87) 

whence the potential function of the second type 1)6000168 
':*n]p=r2 "n-^n"!. ^--7=====^==== !!n(9) * ( (,/fT;i^) ) (88) 

Now the functions Hjj(cp) have orthogonality proper- 
ties similar to the trigonometric functions, as may he 
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proved "by (73) and (74) if it is rememljered that M^(9) 
can "be written in; the form 

Mjj(cp) = (sin cp)^2[ao(cos cp)^"^3~^+ aa(cos cp)^"^^'^ + ...] 

It is ' 

P 2 M^icp) Mjj(cp) Acp d© = 0, if n - "Y = ±2, ± 4, . . . (89) 

The outflow condition reads, according to (88): 
^li M3r («p ) ■* _L Car+i M8r+x(<P) s-O.. . J|<.!p.<. (.90a) 



CO 

r 



which, after multiplication by liza-^ii^) ^ ^ ^''^^ ^--27(9)^9 

respeptively , and integration from cp « to ^9 give 

^sa- i^?a-^;i » 2a-* 1 + 3, ^^r ^srfsa-i = 0 , a = 1,2.,.«. (92a) 

DsV 1-2^.^37 + 2 ^ZT+x I-(sr+i) t-sY 0 V = 1,2,.,. (92h) 
whereby 

. n 

17,6"^/ ^TT M.^^^) M8(cp) Acp dcp (93) 

The infinite equation systems (86) together with (92) 
then enable the determination of the coefficients of the 
potential functions C^^ and D^^ of the second type. 

These coefficients in conjunction with the previously com- 
puted coefficients a^^^ of the functions of the first 
type make the prediction of the aerodynamic q^uantities of 
the given wing possible, as will be Illustrated on a model 
case.. - •• 

Note The calculation is made on the assumption that 
the flow strikes the elliptic disk paraHe! to the minor 
principal axis. The length of the pri ncipal axis in y 
direction was 0, in:- :x direc tion cyi-K^ » The calcula-* 
tion is readily applicable to . the case of elliptic wing 
in flow parallel to the major principal axis when is 
assumed negative, that is, supposing the aerodynamic cen- 
.ter on the y axis. The c;onversion fQvrmulas for the el- 
liptic integrals v;ith negative are: 
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TT TT 

^ d € _ 1 d € 



' 0 



TT TT 



> (94) 



/ / V l^K^ 

-0 -Jo 

The integrals are as far as the factor before the 
integral » the same as for the reciprocal axes ratio. 

The flat ellit)tic disk in straight flow .- let the 
angle of attack he a^, that is / the elliptic area is 
given "by 

z = -X tan » -x 
Then we have, according to (62)f 
w(x,y) * . cLw 

on the disk. According to (83) therefore, the coefficients 

of the potential function of the first type are all 
zero; those of the second type must he computed conform- 
ably to (73, 74, 75a), as exemplified here for n = 3 and 
n = -3. 

For an axis ratio of 

2 



y 1 o . (2c) *t 

Jl - = i = 0.96, that is, A = ^ ^ 6.37 

5 Pell TiVl^K^ 

the complete elliptic integral is: 

TT 

p = /'^ ^ ^ 3.01611; 

Jq VI - 0.96 sin^ € ^ 

3 = J J\ ^ 0.96 sin^€ d€ x 1,05050 
In the case of n « 3 and « .1; =5 ^3 = 0 

E^®(i;) = V'l-r (v^ - Ps"^) s « 1;2 
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and according to (20) it is: 



+ — = 0... that is, pi^ = 0.19792; p3®= 0.97008 



P^-K^ Ps^ Ps^-l 

Conditions (73) and (75a) then read: 
1)3^^'^ / (0.96 sin^e - 0.19792)£.ed£ 

+ ^3 /TT (0.96 sin^ € - 0.97008) cde = 0 



S6 



- Vl - = - ts^^ Vl - K 2(0 . 96 - 0.19792) -1)3 ^/ 1-k2(0. 

-0.97008) 

whence .' ) / ■\ 

= 1.317, 1>3^'"' = 0.3095 

i. e., M3(cp) = 1.627 tt^ sin^ (p - 0.561 

For n = 3 and = ^2 = £3 = 1. we find 



M^*"^( c) = COS € siii e 



Sq.uation (75a) reads : 



.TT TT 



.1 rs >i r sin el) ^ COS € sin c.A e . ^ ^ s ^) sin^ € A € d € 

2 ^ cos € ^ ./ 



that is t 



nence 



1)3"" ^ = 2.654 



(c?) =^-3«&54-co« cp sin cp, J yi^ ^ = - K cos cp (J) 

'. .. ■ ■ ■■ ■ ■ ' ■ 

The determination of the other necessary functions 
proceeds in ■ s imilar manner. The integrals (93) can be com*^ 
puted ly means of the ftmctions JMj^(cp.) ^. "being either ellip- 
tic or reducible to elementary. The coefficients ^y^S 



17ACA Technical Memorandum No. 971 



27 



(92) are herewitJi known. Calculation of the right side 
of (86) yields 

' " r_ ^ ,.^._..H°'o for a=l 



for a >1 



-tto / / . d n' i*sa-i(^)«il= *"^o / ^ ?2a--i(Ti)dTj = \3 

For y = l.S,... the right side of (86) is always 
zero, whence no asymmetrical potential functions occur in 
y. In the effected calculation, the series (81) of the 
potential function was stopped with n = 4; hence equa- 
tions (92) and (86) must he .taken for a = 1; 2. It gives 

2.101 Ci + 1.5217 Cg + 0.527 C4 = 0 (a = 1 )1 

- 0.2410 Ca + 0.5132 C3 + 0.9248 C4 = 0 (a = 2)J 

- i Ci + I Ca + >irC3 * -#ao (a = 1) 1 



f(92a) 



(86) 

20. °i - T¥°3 + f C4 = 0 (a = 2) 



i 



The direct solution gives 

Ci = 0.558 ao, Cgs -0.7785 a© • C3 = -0.342 a© » C4 = -0.0135 a^" 
and the lift, according to (77), at: - 

ii = 0.568 aoXSx^T^ S'ell^**^^ f ^ell 
that is, 

_a 4.55 

a ao 

The moment a'bout the y axis is, according to (78): 

M = -0.7785 x Q.9524 X |.c Vl - £ 7^ p 

^32 ell 

M = -1.98 c yi - k2 £ P^j^j^ 
the center of pressure is at — ^ = -0.435, that is, 

Cv^l - 

at 28.3 percent of the maximum wing chord. 

Incidental to the calculation of th^ induced drag, it 
is emphasized that the lift distribution of the lifting 
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line does not disappear when allowing only for an infinite 
nun*ber of series terms in (81) at the wing tips.(fig- 3), 
and hence must "be included as a sul^stitute lift distribu- 
.tion. In this case the elliptical is most suitable, giv- 
ing for the induced drag the well-known formula 



2 ^ell 
^ n(2c)3^ 



?or the axes ratio 



Vl - ^ =: 0.75, A - 2.55 

the procedure is the same, the quantities "being as- 

certained from (73, 74, 75a), This affords the functions 
Mj^((p) for the integrals ly g, wherewith the coefficients 

of (92) are known. The solution gives 

Ci = 0.3741 a^, Cg = -0.6347 , C 3 -0.2347 , 

C4 = -0.0138 a 

the lift being 

A = 2.9S ao § ^ell ^'^^ 
the pitching moment 



0 



M = -1.397 cjl - ^ Tell 
and the center of pressure at 

X 



= -0.467 or 26.7 percent of maximum wing chord 



For a/i - K^ rr 2 • = - 3 A= 0.637 

it gives 

Ci :r 0,124 a^, Cg =-0.5245 tto , C3 = -0.066 

' C4 « -0.011 ao 

that is , 

d c« 

. : ' . ^—J^ :r 0.99 



UACA : Techaxcsl Memoraudxun Ho, 971 



29 



The center of pressure iff situated at 

— ? = -0.584 a 20.8 percent of the chord 

cvi - k2 * . 

Por comparison the values for the flat circular disk 
are repeated (reference 2) : 



k2 = 1, ^ 0,- ' A = 1.272. = 1.82 

d ao 

center of pressure: — = -O.SIS, i.e., at 24.3 percent of 

c 

■ chord. 

The calculation method used here permi ts even a bound- 
ary transition to the lifting line (k* = vl - = O). 
It affords, v;hen two series terms, are taken into account, 



^ Ca 
d 



= 2 TT c.p. at 28.8 percent of chord, 



The latter result corresponds to an elliptic spanwise 
load distribution with a center of pressure. at 1/4. chord 
in each airfoil section (e.g. of a homogeneous semiellipse) 
Development of all quantities appearing in the equation 
systems vrith respect to affords for small k': 

• . ^ Qa ^ 2 TT . 

d a 



^0 , 6S 
1 + 



128 

d c 

When — - is calcula-ted according to linear wing theory 
d ao :. 

where , as is known , 

1 ' d' Co 2 TT A 

^ ^ TT A d ao A + 2 

there appears a marked discrepancy at small A with re- 
spect to the values computed in accord v;ith the theory of 
the lifting soxrface (fig. 4). • 

On the other hand, the agreement with V/einig's re^ 
suits is good. ( See reference 6. ) - 

Figure 5 shows the center of pressure position plot- 
ted against aspect ration 

The elliptic wing in yaw .- This problem can be treated 
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under the same assumptions as the wing in straight flow* 

The so-called angle of yaw & is defined in figure 6* 

ifow the streamlines are straights in the plane 
z = 0, defined iDy 



(95) 



y « -(x - xr^)tan 3 + \ 
= -X tan p + const -> 

From equation (l) 

(V cos 3 + u) |1 + (-V sin P + y) ^ + w |2L = - i ^ 

OX 3y oz p c^z 

under simplifying assumptions, there is found 

V cos P 1^ - V sin p SI.^V^ |± (96) 
ox cjy C)2 

But "because of (95) it "becomes 

fJi . Sw . SiL tan p\ that is cos P = |il-cos P - ^ sin P 

dx ox dy • dx 3x oy 

Hence equation (96) can te written in the form 

^ cos p = Y ^ (y = - X tan P + const) (96a) 

dx Bz 

The dov/nwash follows from integration along a stream- 
line again assumed as a straight line pariallel to the di- 
rection of flow 

X 



l-r r |i dx (97) 



?. , COS P I bz 



The potential function \if ■ is ag^in assumed as a 
linear combination of the potential function ^-n" 
in the form (8l),. v/hence the sam^e forumlas are obtained 
for lift and moments. 

The dov/nwash function is computed on the "basis of an 
ohl ique-angle system of coordinates in the xy plane, 
given ty the ellipse diameter parallel to the stream direc- 
tion - Ig axis - and the related con jugate diameter 
Tip axis. Posting 
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tan 96 = -^/^"-"^ tan that is , sin 9= = '^^ " sin 3 . 

cos cpp = £^ (98) 

we have the quoted diameters given by 

tan CDo tanKOg + ) 
y = . P X and y = -r respectively (99) 

On the disk, i.e., p = 1, it is according to equa- 
tions (6) and (17); 



I = 1=== = - ^ , cos cp; Ti 5= - = M. sin qp 

c 



_ - - ^ -^ cos o: 5= — 

Ca/1 - yi - k2 

Similarly, we post 



= . cos (cp - cpp) ; r,p = M. sxn (qp - cpg) 

V 1 - k:^ 



that is , 

cos ^ 



yi> 1^ sin P - sin 3 , . cos 3^ ,,^^^ 



In this instance the potential function con- 

formaMe to (32) is utilized. During the respective dif- 
ferentiations and integrations along a streamline, the 
fact that y =s - X tan p + const should be borne in mind. 



//(t) ^ eo X *^ 



- V'e"2"~-'~e^'\/e 



c 



We find that 



(32) 



a 2 d X c - egVe^ - e 



3 



32 



; KACA 'Technical " HenLbr-aixitim. No:. 971 



and ■ • 

iSsi£l= .4^L^ i [_:Zl^lL ^^ ^' ' ^^^^ ^I- tan e 

d X d X c Lv^eg - ©t^eg - © 3 - e^V 63-63 

— S thus can be expressed again "by — 7 ; postxng 

3 z d X 

the result in (97) and integrating over x gives 



P./ ^ Ve, - e^e, - 63 

X ■ 1 .-,r~^^'"(t)dt(lOl) 



^ ^ tan P 

Veg - 6 i^s/es 6 3 3 ©x^y^s ®2 

With equations (51) and (52) is now r94)laced 

again by 

g - 1 r ^n(Y(s)) dY 

The integrand has poles at t = s, where X(t) - T(s) = 0 
and at t - tg» where 

■ -/^(te).j:_e3_' . : ^"^^^P) - ^L. tan P = 0 ' 
*J &z - ei^ 6 2 63 ^/ 63 - e^V^es - e^ 



i, 



e 



A/'Y(tp) - 63 = i Ves. ©3 sin cpp; V'V(tg) - 6.3= -^62 - 63005 cp. 
The denominator at this point is as 

/ ' ' ■ - cos cpe i sin cpg -i 
(t - tp) yv(tp) " 63^ ~ - tan & U 

lIa/ ei 62 ve^ - ©3 J 
The residuum at t « tg is thierefore 
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2niO^(X(t^)) 1 : -S/(tp) 

i cos cpp + i tan cp^ sin <iPp 

= Str i Kcos cpp) <^ (Tip) (102) 

with a Tiew to the fact that, according to (lOO) 
X(tp) = sin cpp — ^ + cos ^ = Tip 



c 



Jl - K 



and according to equation (33) 

E-"'(ta) « i 

P X=- KCOS cpp 

Por th6 residutun at t = 8, the same holds true as 
in the case of straight flow. It affords: 

2Tf i P ip„(|cose +T1 sine) 7===^ ' iMj^°(e)d€ 

/ „ * cos € - -v/l-K^ tan ^ sin € 

5 ___ ^103) 

By replacing a/i - tan P "by tan cpg and posting the 
coordinates |p and T\p "by means of equation (lOO), the 

values (102) and (103) from (lOl) lead to 
^W^«>ap, Tip) = Ea"»(- KCOS cpp)^^(Tip) 



M„"'(0_ 
+cpp) 



(104) 



For the terms of P^, containing powers of |p, the 
denominator cos(c + 9p) in the integral .disappears, af- 
fording complete elliptic integrals of the first and sec- 
ond types. To conpute the others, numerator and denomina- 
tor are expanded with cos(e - cpp)^^ Pt which, with 

cos( e + cpp) cos( € • qpp) p = cos® P - P sin® € 
gives 
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. n _ 

The existent integrals are reducible to complete ellip- 
tic integrals of the first and second types and one integral 
of form ^ 

sin^ c d€ 



cos* 



which, as a complete elliptic integral of the third type, 
is reducible to incomplete first and second types (refer- 
ence 7.) . 

7T • 

2 

(l-K^)sin p A P n sin^C d£ 

cos g / cos^p-A^p sin^c Ac 

= E(P)P(|) - e(|) J(P) - AP tanPF(j^-) (l05) 

The calculation of the downwash function for the po- 
tential function V^^^ is given as a model example. 



E 



Pi(T) = Y q,(X) . I in - 1 



Ei^(- K cos cpg) =5 a/i - cos^ 9^ = ^ ^ 

Then, according to. (104a): 
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1 If ' \ 1 - k2 ^ . . 



TT 



1 ^ / . . , COS (e-^q)Q)A€ 
-A3 / [|p cos(€+ cpg) + Tig sin(c +(Pb)j 1 — ^-T" 



Vl - K 3 ^ > . . 



A«3 



T . 1 1 rs sin9QC0scpft+ sine cose- 

2 Pap 2 P jjjjgs 3 p sinSe 



2 

'^'he evaluation of the integral gives 



v4 



SO)p(l) - s (|) 



The downwash functions of the -potential functions of 
the second tyipe. - According to (68) 



whereljy 



* n m a 



according to (65). The downwash function is 
computed from Wj^® ( ^T^p) , that is, . 

. ii r a F, .(Y) M „ (c) ^^^^ 

2 Ap dT cos(€ + cpg) 

2 

in the same manner. 
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On the disk we have — ^ ' -0, that is, on a 

streamline (tiq = const) the downv/ash function is constant; 
hence w^^ is a function of rip only, and we accordingly 

' ' ' ... I 

put Ift * 0 in ~ ^(Y). To find the summand 

P dT ^ 

Tip^-^Plsin (€ from T^-^P 

P P 

we put T =? ^ and ~7r-~» respectively, whence 

Csin'(c +<Pp)]''''^P =[sin(€ .+ 9p)j'''"^'^[sin(€ + pg)]^'^-^P 

=:[sin(€ + cpp)] (l + sum of cos terms) 

Factor cos(g + 9g) be extracted from the sum, 

thus becoming shorter with respect- to the denominator. The 
new sum, hov/ever, multiplied by [sin(€ +00^)]^^^*^, yields 
only terms which eithei satisfy odd in € or else condi- 
tions (73) and (74), respectively, and accordingly disappear. 
A proportion other than zero is afforded only by 

[Bin(e +9^)]^"^ x 1, which is, however, no longer depend*- 

ent on exponent p. In consequence , all before 
the integral can be combined conformably to 

^ ~— ^n-x^^B^' since the integral for 

dY dT| g 

every p is the same; i.e.. 



7 AP m ^ ^ ' P d Tig 



+ — — ^ ^ I [sine cos cp p + cos € sxn cpg] 
d ^ '-ill 



cos c cos Op + sin € sin 



)c + s in € sm 9a ^ . 

^ M„(€) A € d € (107) 

cos^ P - ^2 g sin^ c 

Hence the downwash functions for the wing in yaw have 
the following form 
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d 



d Tig d Tip 



The second fundamental TJ-rolilein for the case of the 
wing in sidesli-p ,- The procedure is the same as %n the 
case of the wing in horizontal flow, hy putting 

^/T^^vTT^cos pMi!li=I S S g„«» E^^'Cp.) E^"* (v) (108) 

dx c n m ^ 



where 



1 / \ « dz(x,y) 

•| w(x.y) = cos 3 — ' 

and .y,=5 ~3c tan P + const; \|/ to be giirea again by 

= I i ^n" ^n'^ ^l^nK^l^n ^-n 
Equation (96a) then gives ' 

r d p„°(p) 



° ■ a ° 
ci(V^^llJ_ » 



(108a) 



(109) 



(110) 



and with it the coefficients of the potential functions of 
the first type. For the determination of the coefficients 
of the functions of the second type we take (97) in the 
following form: 



n nj 



r 9©H d 



Bz cos 3 



+ S D. 



~>-.oa 



3<|!,n d X 
oz cos 3 



> 



- w(x,y) - S 2. an" / 

V . . n m ./ 



cos P 



(111) 



— ca 



The right-hand side is a function of Ttg only because of 



the determination of the a^ 
rest of the formula then reads: 



from equation (96a) • The 
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Z C k (B) 



d Tip ^ <iTlp 



1 



to 



This equation is integrated from = 0 to ri^ 

then - multiplied by ^sa-i^^g^ ^^sVCnp) again 

from 'r\^ = -1 to +1» This affords an the basis of the 
orthogonality characteristics of the spherical functions 



and 



2 (C; 



r«o 



sr+i^sr+i 



0)+B 



2r+i ^sr+i 



a(2a-l)-r(2r+l) 
2 



1 



4a -1 



1 
T 



I i 

I 



J . 

— X 0 



w(n 'p) dTi 'g?2a_i (Tip ) dTi^ 



ex ~ X \ 2 f • m • 



(113a) 



2 
r=x 



(Car '^BT (&) + 32r Isr O)) r 



y(27 - 1) - r(2r - l) 
1 



3^ 
V 



47+1 



(113b) 



The outflow condition ^- The formulation of the outflow 
condition is prefaced by the following note. The downwash 
functions (107) on the wing in sideslip have the quality of 
becoming infinite by T^g = ±1 like ^ ,(Tia) and 

, respectively. But Tj^. =s ±1 are marginal 

d Tip ' p.. 

points of the elliptic disk in which parallels to the flow 
direction touch the ellipse (fig, 6). Accordingly, it is 
to be assumed that the sorcalled "vortex tails " in these 
•points leave the disk parallel to the direction of the 
stream, Eence it is postulated that no flow around the 
trailing edge occurs between the points* Tjg = ±1, In 
other words ; 
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This equation, multiplied /by 
Isa-i^^^ ^ ^ and M^gV^^p) respectively 

and integrated from 5" + ^Pp ^ + aI^q gives the f ol- 

2 P 2 P 

lowing systems of equations: 

P tO) J_ V P T^P) . r-. T\ T^P^ r\ 

CO ~ ly2f*** ^ 

wherehy 

8 = / HvCqp) MgCqp) Acp d(p (114a) 

These two equations ((113) and (114)) make it possihle 
to determine the coefficients and J>^. Together with 

the previously defined a^^"^ 's the aerodynamic quantities 

of an elliptic wing in yaw can he computed. As to the 
equation systems themselves, they form a coupling of the 
systems (82) and (92) for the case of straight flow, as is 
readily apparent from the similarity of the corresponding 
coefficients and which thus affords a first simple mathe- 
matical check. 

The flat elli-ptic wing in yaw .-- The calculations are 
carried out for the axes ratio 



5 

and the angles of yaw 

e = 15^ and 3 30^ 

Again only the potential fvmctions of the second type 
conformable to (108) and (llO) are required. To simplify 
the voluminous paperwork only the functions up to the de- 
gree n = 3 are taken, i.e., 
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^'Xf ^39 ^^B9 • (^-1 ^O^.S ^Ot exist) 

The five unknown coefficients of these potential func- 
tions require five equations^ Two are taken from the con- 
dition i w(T|g) = - cos p tan on the disk, equation 
(113a) for a = 1, and (I13b) for Y = 1. For the other 
three the outflow condition (equation (114)) with a- 1*2 
and 'Y r 1) is used. The Mnv9) are those previously de- 
fined in the case of straight flow. The integrals 

I^P^ ostein now only the limits f + and + cpg, 

while 1^ y remain unchanged on account of the periodicity 
of the integrand* The values of the incomplete elliptic 
integrals of the first and second types necessary for the 
determination of the coefficients k^(p), ^n^^^* etc;; 
v/ere taken from Legendre's tables (reference 8) 



i • e • , 



For 


^ 


= 15° 


it affords 






Cx = 


0.5204 






Cs = 


-0.7194 ttjj = 0.0144 








-0.3343 ttQ D3 = -0,0065 






A = 


4.16 1 S-^^i 






M = • 


-1.83 ctQ c yi - k2 1 V» Pg^j^ 






L = 


0.0384 ao C | S'ell 


Pqr 


.p. 


= S0° 


it is; 






Ci = 


0.407 ttQ 






Ca = 


-0,562 ttQ D3 = 0,0277 



i.e.. 



C3 = -0.265 = -0,0124 a. 



A = 3.26 tto £ Pgj^i 

M = -1.43 c -I Pgj^j^ 

L » 0.074 tto c I V2 P^ii 
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The new additive moment L is positive according to 
the above , calculations , which is synonymous with the fact 
that the leading wing half receives greater lift,- The 
coordinates of the centers of pressure are: 

p = I5Q; — ^ -0.440; - = 0.00925 = 6^020 

c^l - ® 

P = 30^1 ^ ,.. = -0.439; I = 0.0227 (cp = 13^510 

For comparison we repeat the values obtained at p = 0^ 
when the expansion is stopped with n s= 3. The bracketed 
terms contain the change in percent vflth r'espect to the 
quantities computed with the four expansion terms, and from 
which inferences can be made regarding the convergence. 

Ci = 0.563 ao (-0.9 percent) 

C2 = -0.776 (-0.2 percent) D2 = 0 

C3 = -0.365 tto (+7.0 percent) = 0 



i.e., 



A = 4.50 tto £ V® 



2 ell 



M = -1.98 ao cyr^-K^ I V2 V^^^ 

Center of pressure: — ^ = - 0.438 (+0.7 percent) 

c/l - 

The results are correlated in figures 7 and 8. 

For great A an approximate formula for the rolling 
moment in relation to angle of yaw ^ and axes ratio 



K • = V 1 «. is again expedient: 



lntan(j+i) 



In i ^ 

j£L,ii.^3 K' sin 2^ ° '^'"2" iM-J ^ (115) 

io-o „ 243 

1 + Kl 

cos ^ 138 

The formula indicates that the moment on transition 
to the lifting line (k' = O) disappears. But if the 
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moment with the half wing ^ chord: ins tead. .of half the span 
is made nondimens ional , thus voiding the factor in 
(115) the moment coefficient "becomes' logarithmically in- 
finite on limiting transition. In the extreme case the 
lift decreases with cos^ for the reason that the 

flow velocity in x direction- is V ;c;as ' 

Horner's results on v;ings of different plan forms in 
sideslip are in very close af:r^ement with the values given 
here. The assumption that the rolling moment is, aside 
from the angle of yav;, largely dependent upon the aspect 
ratio rather than the chord distribution appears, therefore 
justified. 

. -J 

Translatioii "by J, Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs, 1,2, 3 




Fig.l Base Ellipsoid 
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Computei lift distribution of the 
lifting line. 

Elliptic lift distriliuti,on with 
equal total lift. 

Fig, 3 
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Fig-B. 4,5 
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Figure 5. 



iJACA Technical MeBorandum No. 971 



Fiss.6,7,8 




Figure 6. 
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Figure 7 plotted against angle of yaw p 

for aspect ratio y\. =6.37 
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Figure 8. ^ plotted affainst angle of yaw 3 

da. 

for aspect ratio/V = 5.37, 



